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j5 . Abstract. Recently, Nattermann and Pokrovsky [PRL 100, 060402 (2008)] have 

^ ' proposed a scaling approach for studying Bose-Einstein condensates in strongly 

"^-J '■ disordered traps. In this paper we implement their scaling argument in the framework 

, of the variational method for solving the time dependent Gross-Pitaevskii equation. 

We consider atomic gases with both short range s-wave interaction and long range 
anisotropic dipolar interaction. The theory is addressed to the regime of strong disorder 
' ^ I and weak interactions where the physics is dominated by the collective pinning due 

Q ' to the disorder. The phenomenon of condensate fragmentation in dipolar gases is also 

O . analyzed. 
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Introduction 

Motivated by tlie recent progresses in ultracold alkali atoms in disordered traps [21 |3l H] , 
Nattermann and Pokrovsky proposed [1] a semiquantitative approach for Bose-Einstein 
condensates in strong random potentials where standard perturbative methods are 
expected to fail. Their analysis is based on the evaluation of the mean-field energy due 
to the fluctuations of the random potential according to a scaling argument introduced 
by Larkin [5] for studying the effect of defects in flux line lattices. The scaling 
approximation requires that the correlation length of the disorder is much shorter 
than the correlation length of the liquid. In this paper, we implement the method 
of Nattermann and Pokrovsky in the framework of the time dependent variational 
method of Perez-Garcia et al. [6]. As a result, we obtain a simple set of equations, 
which describe the equilibrium and the low-energy dynamics at zero temperature of a 
Bose-Einstein condensate in a strong random potential. Beside the oscillator and the 
scattering length, a new quantity enters into the problem, namely the Larkin length £ 
associated with the collective pinning of the condensate. 

The time-dependent variational method represents a very good approximation of 
the Gross-Pitaevskii equation [TJ Ej. Moreover, it can interpolate quite successfully 
from the low density regime to the strong coupling Thomas-Fermi gas. The present 
extension to disordered traps, however, presents some restrictions. Since the fluctuating 
center of attraction of the disorder may not coincide with the center of the harmonic 
trap, the model is fully consistent only when either the harmonic trap or the disorder is 
responsible for the localization of the condensate. Moreover, the Gaussian variational 
ansatz out-rules from the outset the solutions describing the multi-fragmented state 
of the condensate [1]. Despite these limitations, novel quantitative predictions can be 
obtained with regard to the disorder dominated regime. In this limit the condensate 
becomes non-superfiuid and is characterized by a generalized correlation length [9] larger 
than the size of the cloud. 

The paper is organized as follows. In Section 11.11 we introduce the general 
formulation of the problem. We determine the size of the ground state at equilibrium and 
we calculate the expressions for the frequencies of the low-lying excitations. In Section 
11.21 we discuss in detail the 3D gas in presence of a strong delta-correlated disorder. 
We show that, reducing the atomic scattering length by means of Feshbach resonance 
techniques, the predictions of the theory for the single connected fragment might be 
tested experimentally. In Section 11.31 low-dimensional gases are considered. Section 11.41 
contains some considerations on the applicability of the method in the case of disorder 
with finite correlation length. In Section 12.11 we introduce a long-range anisotropic 
dipolar interaction and we study the stability of the single connected fragment. The 
nature of fragmentation in presence of dipolar interaction is analyzed in Section 12.21 
The static on-average characteristic properties of the fragments are determined at the 
level of the Imry-Ma level of approximation [10] . 
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1. Time dependent variational method with disorder 
1.1. Variational equations 

We begin considering the Hamiltonian of a gas of trapped bosons interacting through 
s-wave scattering 

where a is the s-wave scattering length and the harmonic trapping potential is Virap(x) = 
^i=x y z''^^i^'i f^- The disorder potential f/(x) is chosen to be Gaussian distributed 
characterized by the average values (U) = and (f/(x)f/(x')) = (K^/fo"^) (x — x'), 
where b denotes the correlation length of the disorder. We will assume that -^'^(x) is a 
smeared out (5-function -K'b(x) = e~''^/^^/(27r)'^/^. 

The variational calculation follows the same outline of the clean case. Therefore, in 
what follows, we refer the reader to [B] for the details of the derivation. We consider the 
normalized variational wave-function in the case of a 3-D fully anisotropic configuration 



hiBiit) 



where we have defined the geometric average R = [Ra:{t) Ry{t) Rz{t)]^^^ and the variables 
Ri and Bi are time dependent variational parameters. The i?i's are related to size of 
the system while the imaginary width iBi of the Gaussian ansatz is necessary to include 
the dynamics into the variational procedure. Without this latter the minimum principle 
would just lead to the condition of static equilibrium. 

Physically, the true ground state of the stationary Gross-Pitaevskii equation in the 
presence of interactions and disorder is expected to deviate strongly from the ground 
state of the noninteracting system. Nevertheless, for a clean system, it has been shown 
that the dynamical trial wave-function ([2]) is a very good approximation of the Gross- 
Pitaevskii at finite [6l [71 E]. Inserting the trial- wave function in the Lagrangian 
relative to the Hamiltonian of ([1]) we find 

K iV 

where Lq is the well known contribution in absence of disorder [6]. The new term 
represents the on-average contribution of the disorder fluctuations according to the 
scaling argument of Nattermann and Pokrovsky [H Ej. It introduces a new relevant 
length scale, namely the Larkin length C = (2^/^7r'^/^/i^/m^K^), associated with the 
pinning energy due to the disorder [1^. The length scale £ was introduced by 
Larkin [5] in connection with the onset of collective pinning of vortex lines in type- 
II superconductors. In the case of a trapped gas, the analogous of the (delocalizing) 
elastic energy relative to the lattice distortion corresponds to the kinetic energy of the 
atoms. Since the scaling argument requires to have both C and b much smaller than 
the generalized healing length ^heai [9], z.e. £ <^ ^heai and b <^ ^heai, the validity of ([3]) 
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is in general restricted to the quantum limit b < C. The classical limit b ^ C, when 
many levels occupy the typical well, will be briefly discussed separately. 

In what follows it is convenient to rescale the quantities in unities of the harmonic 
oscillator such as rj = Ri/i and (b/ri) = {b/i){i/Ri) where i = ^Jh/muJ is the harmonic 
oscillator length. The anisotropy of the external trapping is taken into account by setting 
cui = AjCU. Without loss of generality we will consider a trap with cylindric symmetry 
with anisotropy factors A^. = = 1 and A^ = A. In absence of disorder this implies 
that the angular momentum along the z axis is conserved and we can label the modes 
by the azimuthal angular quantum numbers m. 

Next, we derive the Euler-Lagrange equations for the variables Ri{t) and Bi{t) 
relative to the Lagrangian ([3]). Then, eliminating the -Bj's variables, we obtain a closed 
system of differential equations for the i?j's. These equations can be viewed as describing 
the motion of a point-like particle in an effective potential. The position of equilibrium 
at rest of the particle is determined by the minima of the effective potential. Using 
rescaled unities r^^ = and tqz, this leads ultimately to the conditions 
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where we have deflned the functions fi {tq. 
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and gi (roi) = 



b/ro, 



fi/2, and we have introduced the Thomas-Fermi parameter a = \/2/7r (Na/i) 

together with the disorder strength parameter 7 = ^JTJL. In these two latter equations, 
the l.h.s. describes the conflnement due to the trap and to the disorder. In the r.h.s. this 
effect is counterbalanced by the kinetic energy and the repulsive mean-fleld interactions. 
The limit £ — ?■ 00 corresponds to the theory for a clean system [6]. 

The low-lying excitations of the system can be calculated by means of a harmonic 
expansion rj(t) = -|- 5rj(t) around the equilibrium position (jl]) and ([5]). Making a 
dynamical ansatz for the 6ri{t), the calculus of the frequencies for the small oscillations 
is reduced to the calculation of the eigenvalues of a 3 x 3 real symmetric matrix Aj^ 
whose elements are 
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We denote the eigenvalues as cUabc- In absence of disorder 
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corresponds to the 
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Figure 1. (a) Frequency of the breathing mode Wb (red sohd hne) and of the 
quadrupole e-mode (blue sohd line) as a function of the particle number. The upper 
horizontal green dashed lines describe the analytic results of ([7])- (|8]) while the lower 
dashed magenta correspond to the Thomas- Fermi limit ujh = \/5w and ujc = \/2uj. We 
have considered a sample of ^''Rb atoms confined in an isotropic trap with frequency 
w = 27r50 Hz and with scattering length a = SOag, where oq is the atomic Bohr 
radius. We have considered delta correlated disorder with C ~ 200a which amounts 
to £/C ~ 2.9. (b) Same frequencies as in (a) and the ratio LUf,/iOc (magenta line), for a 
very shallow trap with uj = 27rl0 Hz and for a very small value of the scattering length 
a = lOao- The Larkin length is £ = 2000a which means £/C « 3.2. 

frequencies of the quadrupole modes with quantum numbers n = and / = 2, while 
is the frequency of the monopole mode with n = 1, I = 0. 

It is important to remember that the scaling argument leading to (|3]) applies under 
the strong disorder condition i ^ C Since the center of the attractive domain formed 
by the disorder may not coincide with the center of the harmonic trap, the variational 
ground state is exact only when the system is trapped either by the static fluctuations of 
the random potential or by the harmonic trap. In principle, one could try to extend the 
ansatz ([2]) including an additional variational parameter which describes the center of 
the condensate. However, the disorder term in ([3]) would not bring any new contributions 
to the equation of motion of the new variable with respect to that of the clean model ^ . 
This shows that the interplay between the two different confining mechanisms cannot be 
incorporated in the theory, which is a direct consequence of the on-average nature of the 
Larkin energy considered in ([3]). Therefore, the crossover region has to be considered at 
best as an extrapolation. 

Moreover, our ansatz for the wave-function cannot describe a multi-fragmented 
state, which is expected to occur in a wide range of parameters [H [12]. In this case 
our solution may occur as a metastable state [1]. Due to these limitations, in the next 
sections we will focus our analysis mainly to the single fragment non-superfluid regime 
where the cloud is trapped by the disorder and the physics is dominated by the Larkin 
length. 
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1.2. Strong disorder with zero correlation length 

In the limit of very short correlated disorder b < C the corrections due to finite 
correlation length are very small. Therefore, as first approximation, we can consider 
a delta correlated disorder, and put b = 0, fi = 1 and gi = in ([3])- ([6]). Moreover, in 
this limit and for low energy oscillations, the wave-length of the modes is much larger 
than the distance over which the disorder varies. Therefore, the oscillations can be 
considered self-averaging and we expect the values of the frequencies to be independent 
of the specific realization. 

At low densities and strong disorder, when both the radial and the axial oscillator 
lengths of the traps are larger than the Larkin length, the gas is confined mainly by the 
random potential and the cloud is spherically symmetric. The equilibrium conditions 
in (jlj) and ([5]), which determines the size of the system, can be approximated by 
~ 1 + a/ro. 

In the limit N <^ {^/o) the interactions can be neglected, the size of the cloud is i? ~ £ 
[13] and we find 

UJa = U}k= CQUj{i/Cf (7) 

uJ, = {co|2)uJ{^/C)\ (8) 

where the constant cq depends on the normalization of the variational trial wave- 
function we have chosen. Nevertheless, although the exact physical value of the constant 
Co cannot be rigorously determined by our approach, the ratio uj^/oJc = 2 of the 
quadrupole and the monopole does not depend on it. For larger number of particles 
C<^Na<^ {i/Cf^i the size of the system is determined by the competition between 
disorder and interaction and we find 
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In this interval the generalized healing length ^ reaches at some point the size of 
the system which is of the order of {Na^^^ C}^^ . Moreover, the two curves exhibit 
an "avoided crossing-like" feature. For even larger number of particles the external 
harmonic potential dominates over the disorder and the frequencies approach the well- 
known Thomas- Fermi [15]. The crossover between the different regimes of the values 
of the quadrupole and monopole modes is illustrated in figure [1^ for some typical 
experimental parameters. 

At iV > [ilCf'\ila) the crossover to the Thomas-Fermi regime of the clean 
theory occurs [1]. Deep into the Thomas- Fermi regime the system is superfiuid with 
Cheai ^ C.- The critical number of particles where transition to superfiuid occurs 
can be estimated by equating the Larkin length with the superfiuid healing length 
obtained from eqs. (!l])-(l5l) and the definition given in [9]. This leads to A^c ~ jaC? in 
agreement with the theory of [131 [21] . In the Thomas- Fermi limit the disorder can be 
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Figure 2. (Colour online) (a) Example of the oscillations frequencies as function of 
the number of particles for negative scattering length. The upper solid (green) line 
describes the a, 6— modes while the lower (red) shows the c— mode. The dashed (yellow) 
lines indicates the asymptotic solutions (I?])-®, (b) Frequency of "quadrupole" mode c 
(red solid line). The upper dot-dashed light blue line corresponds to the analytic result 
of (fT2)) . The solid dark-blue line indicates the frequency a/3 a;|| characteristic of the 
ID mean- field theory, while the solid light-blue represents the A ^> limit of the 3D 
Thomas-Fermi frequency \/5/2 W|| . The clean theory is described by the dashed green 
line which at low density approaches the free h.o. value 2a;|| (magenta dot-dashed). 
The inset shows the details of the crossover from the ID mean-field regime to the 
Thomas-Fermi condensate. 



treated perturbatively and a small negative linear shift in ^hca.i/'^ would be expected [T6] . 
However, in figure [1^ a positive shift of the frequency of the breathing mode is found. 
This is not surprising since the present theory has to be considered as exact only in the 
opposite limit ^heai ^ where disorder is dominating. 

The frequencies plotted in figure refer to a single connected condensate. However, 
above N > C/a the single connected condensate is a higher energy metastable state, 
since the ground state is expected to undergo fragmentation [1]. In the fragmented state, 
we do not expect sharply defined frequencies. More likely, they should be distributed in 
some interval. Note that the possibility to observe the single connected condensate, in 
the regime where fragmentation is expected, could arise by means of a sudden decrease 
of the oscillator frequency. However, it is difficult to make further estimates about the 
different timescales involved in the experiment. In order to prevent fragmentation, we 
would like to push the single fragment solution at larger particle numbers. For a given 
disorder (£) and trap configuration (£) such that strong disorder condition £ ^ £ is 
satisfied, this can be achieved by maximizing the ratio C/a. Experimentally, this can 
be obtained by tuning the scattering length close to zero [zero crossing) by means of a 
Feshbach resonance as proposed in Fig. [T]d. 

In the case of negative scattering length there is only a metastable state of finite 
radius. This state becomes unstable at a critical particle number Nc- In Fig. [2ti, we 
plot the frequencies of the three modes in the case of an attractive interaction for an 
isotropic trap with u = 27r25 Hz, a = — lOOao and £ = 200a. The divergence signals 
the instability of the gas when the number of particles reaches the critical value N^.. 
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1.3. Anisotropic traps and lower dimensions 

We have seen that, at low densities, when the physics is determined only by the balance 
between the random potential and the kinetic energy, the frequencies of the oscillations 
do not depend on the external harmonic potential. Therefore, the anisotropy of the 
trap becomes irrelevant. Increasing the number of particles, the mean-field interaction 
becomes important and the harmonic trap affects indirectly the results through the 
Thomas-Fermi parameter a. At even larger particle number the system enters the 
Thomas-Fermi regime where the harmonic trap dominates over the disorder and the 
anisotropy plays its usual role [15]. If we consider, for example, an elongated cigar 
along the axial 2;-direction, with Ux,y = a;_L ^ uj\\ = Uz, this description holds as long as 
we have £ <C £_l -C . However, when <^ C i\\ the system can pass through the 
different regimes in the axial directions while its ground state in the radial direction is 
frozen to that of the harmonic oscillator of the radial harmonic confinement. 

Herewith, we limit our discussion to the single connected solution, neglecting 
the fragmented low-dimensional condensate where it occurs [1]. Using the definitions 
o-iD = ^i/o and Cid = (£^1)^^'^ have that for small N <C aic/Cin the axial radius 
is determined by the equation 7r^Q^ = 1. This leads to the results 

UJa = UJb = 2UJ^ (11) 

io,ocuj^{iJCf^^ . (12) 

At larger such am/ Cm ^ N <^ am/ CiD{i\\/ CmY^^ the axial radius is determined 
in the leading order by the equation 7r^Q^ = ar^o and we have in first approximation 

uja = Ub = (13) 
uj, oc uj^ {IJCf {IjNaf^ . (14) 

At even larger values the Thomas-Fermi regime is realized in the axial direction while 
in the radial direction the ground state wave-function can still be determined by the 
free harmonic oscillator. Therefore, we have that the frequency of the breathing mode 
approaches the value Uc = a/3w|| characteristic of the so-called ID mean-field regime 
p!7] . In order to appear, this regime requires the two conditions NaX/i± <^ 1 and 
{Na/y/Xi±y/^ > 1 [18]. Finally, for NaX/£j_ > 1 the gas enters the full Thomas-Fermi 
regime and the "quadrupole" mode approaches the strongly anisotropic 3D Thomas- 
Fermi result [15] ooc = \/5j2uj\\ while Ua = V2oo± and oob = 2uj±. The full crossover 
is shown in Fig. [2]d for the lowest axial mode c, in a ^^Rb condensate confined in a 
strongly anisotropic cigar-shaped trap with u± = 27rl50Hz, uj\\ = 0.01u±, a = 150ao, 
and £ = 300a. For these parameters, i_L/C ~ 0.36, and i\\/C ~ 3.7. Moreover, the 
gas cannot enter the strong interacting Tonks regime since it would require to fulfill 
the two conditions ij_/a < 10^ and A^A < 1 simultaneously [18]. In our case, we have 
\f\i^/a ^11. Nevertheless, the Tonks gas remains beyond the reach of our mean-field 
approach. 

The treatment of a 2D disc-like geometry trap when A oo follows essentially 
the same outline. In this case we define io^^y = ^ = w^. Similar to the ID 
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configuration, a non trivial interplay between disorder and low-dimensionality arises 
only when i±_ <^ C <^ Except that for very large N, the ground state in z-radial 
direction coincides with the lowest eigenstate of the harmonic trap, namely = 
The 2D analog of (fT2!) and (fT4|) are respectively Uc oc u±l±/C in the small limit 
and Uc oc {l±/ C){1_l/ NaY^"^ in the moderate interacting regime. At larger values of N, 
for sufficiently strong anisotropies, the size of the cloud in the x, y plane is determined 
by the competition between harmonic trap and interactions while in the z-direction the 
ground state can still be that of the free harmonic oscillator. This is the 2D analog of 
the ID mean- field regime for which Uc = 2u\\ [19]. Finally when N oo the A — oo 
limit of the 3D Thomas- Fermi is recovered where ujc = ^/wjs uj\\ [T5] . 

1.4- Finite correlation length 

The delta correlated disorder approximation discussed so far, is appropriate for disorder 
with very short correlation length, i.e. when b < C. In this case, the correlation length 
of the disorder is generally much shorter than the correlation length of the system, which 
represents a necessary condition for applying the Larkin scaling argument. The small 
corrections due to finite b can be taken into account using the full equations (jl])-©. 
Differently, when b ^ C, the use of the expression for disorder energy in ([3]) is somehow 
questionable. Moreover, also self-averaging is expected to be much less efficient than in 
the limit delta correlated disorder. Therefore, we conclude that the description based 
on (jll)-(l6]) holds as far as the limiting case of delta correlated disorder is approached. 

Nevertheless, it is interesting to make the following remark. The classical limit 
b ^ C has been recently investigated in [lH [20] . Near one of its typical minima, the 
random potential can be approximated by a harmonic potential well of depth Uq = n/b^/"^ 
and width b. The bound state is located very close to the minima of the random potential 
and the number of levels in a well is large. From Uq = mubb'^/2 we obtain the oscillator 
length of the typical fiuctuating well as 4 = 6 (C/b)^^^ [211 120]. Interestingly, this length 
can be obtained from (j4]) in the limit of small number of particle where the harmonic 
trap and the mean-field energy can be neglected. At larger particle number the size 
of the fragment is determined by the competition of the mean-field interaction and the 
disorder energy. This yields R ~ {NaC^/H'^/^y/^, which agrees with the size of the 
fragments in the multi-domain state found in [21] by different methods. This analysis 
reveals that the the Gaussian ansatz ([2]) is able to reproduce the ground-state localized 
state also for disorder with long-range correlated potential. 

2. Dipolar gas 

2.1. Stability of the gas 

The time dependent variational approach can be extended to the presence of the long 
range anisotropic dipolar interaction [22]. If we assume the dipoles aligned along the 
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Figure 3. (Colour onliiic)(a) The size of the single fragment of ^^Cr atoms along the 
radial (green) and the axial (magenta) directions as a function of the "relative" dipolar 
strength. The parameters are N = 1000, a = IOoq, A = 1, = 27rl0 Hz, C ^ 2000a, 
and i/£ ~ 3.2. The inset shows their ratio ro^/roz- (b) Dependence of critical value 
of the dipole interaction (odip/a) from the particle number for C ^ 2000a (red) and 
for C ~ 1000a (green). The other parameters are taken as in figure [3lja). The dashed 
line at = 150 is just a guide for the eye. 



^-direction, the interaction between two polarized dipoles can be written as 

■2 f¥Y,o{9) 



/io/i' 



4 



4n -y 5 r3 ' ^^^^ 
where 6 is the angle between r and the direction along which the dipoles are pointing. 
The constant fiQ is the magnetic constant and fi is the magnetic moment of the 
atoms. The strength of the dipole interaction relative to the short range potential 
will be expressed through the dimensionless quantity e = ^/^lJb {^Q^'^m/ Air'^h^a). 
Alternatively, the strength of the two potentials can be compared by defining a 
characteristic length scale associate to the dipole-dipole interaction. From the 
uncertainty principle, we can define Odip as the distance at which the dipolar potential 
energy ( 1T5|) equals the kinetic energy. This yields 



ip 



(16) 
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where we have chosen the prefactor such that e = adip/a [23] . 

In presence of disorder and dipolar interaction, the equilibrium conditions of (jHIS]) 
become 
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with |22] J^(0 = [v^/6(l- 
[V^/3{l-ef] [-2^^ + 10^2 + 1 



ey][-Ae-7e + 2+9enm, and 

- 9^^n (0], and n (0 = Arctanhv^r 
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Analogously to the clean gas, these equations have at most only one stable solution. Here 
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Figure 4. (Colour online) (a) Energy of a ^^Cr condensate in a strong random potential 
in absence of dipolar interaction (red) and for Odip = 0.40a (blue) as function of N. 
The other parameters are a = IOOoq, A = 1, a; = 27r50 Hz, £ ~ 30a, and £/C ^ 9.5. 
The inset shows the correspondent energy per particle, (b) Typical deformation of the 
fragments as function of the relative dipolar strength for two different values of the 
disorder C ^ 50a (green), C ^ 30a (red). The other parameters have been chosen as in 
figureHK. ^^Cr atoms, a = lOOao, A = 1, w = 27r50 Hz, upper C ^ 50a, lower C ^ 30a. 

we are interested in the regime of strong disorder and moderate interactions where the 
gas is confined mainly by the random potential and is constituted by a single connected 
fragment. In this limit the anisotropy of the harmonic trap does not play any role and, 
in absence of dipolar interaction, the cloud would have spherical symmetry. However, 
turning on the dipolar potential the fragment tends to become more prolate along the 
direction of the dipoles orientation. The situation is illustrated in figure [HK, where we 
have considered a very short scattering length in order to shift towards large N the onset 
of multiple fragmentation which occurs above > C/a. Moreover, such a small value 
of the scattering length, permits to access the strong dipolar regime Odip/a > 1 even for 
gases made of atoms with small magnetic moments. This regime has been recently real- 
ized in experiments with ^^Cr condensates near a Feshbach resonance [251 [26]. For any 
fixed N, the fragment collapses at some critical value of Cdip/a. In figure [3b we plot this 
threshold as function of the number of particles at fixed a for two different values of the 
Larkin length. For given a and at shorter Larkin length the condensate is less stable 
against the attractive dipole interaction. Note that a shorter Larkin length amounts to 
a stronger center of attraction due to the fiuctuations of the random potential and thus 
to a smaller volume occupied by the fragment. 

2.2. Fragmented state 

So far we have not discussed the possibility of a multi-fragmented state since the 
variational method is based on a single connected ground state. However, in [1] it has 
been shown that for a large region of the parameters space, strong disorder favours 
fragmentation. It is therefore interesting to analyze the situation when a dipolar 
interaction is present as well. In the spirit of [l], we evaluate the Hamiltonian operator 
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on the ground state ([2]). Assuming cylindrical symmetry and rescaling energies in unities 
of the harmonic trap energy hw, the energy per particle of the system can be written as 

in terms of the variational parameters r = r^ and cr = r/r^. The deformation function 
hiO = v^Stt [1 + 2^'^ — {^)] /3 (1 — ^'^) is continuous, monotonous and positive for 
^ < 1 and negative for ^ > 1. 

In absence of dipolar interaction and for a strong given disorder potential the total 
energy E{N) = N£ has a minimum Ei{Ni) at negative energy as shown in figure H^. 
According to the Imry-Ma criterium [l] , at larger than A^i it is energetically favourable 
to split the gas into fragments of energy Ei. Eventually, at even larger A^, the oscillator 
energy per particle of the fragmented state reaches the energy per particle from the 
disorder and the crossover to the Thomas- Fermi regime occurs [1]. 

The effect of the dipolar interaction on the phenomenon of the fragmentation can 
be seen as follows. The attractive dipole interaction shrinks the total condensate volume 
causing an increase of the kinetic energy which tends to delocalize the condensate against 
the collective pinning induced by the disorder. This is consistent with the numerical 
results shown in figure H^. At fixed disorder, finite values of the dipolar interaction shift 
upward the minimum of the energy E = NS until it becomes positive. Then, strictly 
speaking, the possibility of having fragmentation is ruled out at least at the level of the 
Imry-Ma argument. Nevertheless, at sufficiently small dipolar interaction fragmentation 
can occur. In the presence of anisotropic interactions the Imry-Ma [10] argument does 
not imply that domains have to be spherical. The typical fragment will be preferentially 
cigar-shaped with the long axis parallel to the dipoles |27]. This means that, while in a 
single realization of the disorder the various fragments should appear with random shape, 
they will have on average a characteristic deformation. The dependence on the dipolar 
interaction of the typical deformation can be estimated evaluating the equilibrium radii 
ad rQz in the minimum of the energy curve E{N). The result obtained is shown 
in figure Hb for two different realization of the disorder. The curves stop where the 
Imry-Ma argument for the existence of fragmentation becomes invalid. Note that, for 
the chosen parameters, this boundary is smaller than the critical dipolar strength at 
which the cloud would collapse. 

Conclusions 

We have proposed a hybrid approach, originated from the variational method for the 
Gross-Pitaevskii theory and the Larkin-Imry-Ma scaling argument of [1], aimed to study 
zero temperature Bose-Einstein condensates in strongly disordered traps. Similarly as 
for other disordered systems, the Larkin-Imry-Ma scaling analysis, despite its intrinsic 
limitations, allows to make semi-quantitative predictions in situations where no mean- 
field description seems possible |28] . The theory addresses the problem of the condensate 
in the limit of strong disorder and moderate interaction. There, the dynamics of the 
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condensate is dominated by a new length scale, namely the Larkin length, associated 
to the collective pinning of the condensate. In the case of isotropic interaction, we 
have investigated the stability and the low-lying excitations of the system at different 
dimensionalities when the ground state consists of a single connected condensate. We 
have suggested that, for given disorder strength and trapping parameters, this regime 
could be observed experimentally by using scattering length zero crossing [29]. In the 
case of dipolar interactions, our theory offers in addition the possibility to investigate the 
effects of the anisotropic interaction on the phenomenon of fragmentation. Predictions 
about the typical deformation ratio of the fragments can be made. We conclude that, 
the realization of types of disorder with correlation length smaller than all other length 
scales, could allow to test the theory in the laboratory. 
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